Statistical mechanics of an ideal Bose gas in a confined geometry 
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Abstract 

We study the behaviour of an ideal non-relativistic Bose gas in a three-dimensional space where 
one of the dimensions is compactified to form a circle. In this case there is no phase transition 
like that for the case of an infinite volume, nevertheless Bose-Einstein condensation signified by a 
sudden buildup of particles in the ground state can occur. We use the grand canonical ensemble to 
study this problem. In particular, the specific heat is evaluated numerically, as well as analytically 
in certain limits. We show analytically how the familiar result for the specific heat is recovered 
as we let the size of the circle become large so that the infinite volume limit is approached. We 
also examine in detail the behaviour of the chemical potential and establish the precise manner in 
which it approaches zero as the volume becomes large. 
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I. INTRODUCTION 



The study of the idea 
mechanics. (See Refs. 



Bose gas is a standard topic in any treatment of quantum statistical 



HQS 



for example.) The usual assumption made is that the bulk, 
or thermodynamic, limit is taken where the total number of particles and volume of the 
system are taken to infinity with the particle density fixed. In such a limit, some of the 
thermodynamic quantities that describe the system, such as the specific heat at constant 
volume, are not smooth functions of the temperature indicating that a phase transition 
occurs. This phase transition is usually interpreted as Bose-Einstein condensation of the 
gas. 

If we do not take the bulk limit, by keeping the volume of the system finite, then ther- 
modynamic quantities like the specific heat remain smooth functions of temperature. (This 
was emphasized by van Hove for instance.) This smooth behaviour has been borne out in 
a number of different studies (see the review of Pathria j^] and references therein), and will 
be shown below. Nevertheless, it seems intuitively plausible that as we let the volume of the 
system become very large, but still finite, the difference between the non-smooth behaviour 
in the bulk system and the smooth behaviour in the finite system should become less. The 
principal aim of this paper is to study in some detail using analytical methods how the bulk 
limit is approached. 

In order to make the calculations as simple as possible we will limit our attention to a 
three-dimensional space with only one compact dimension. The compact dimension will be 
chosen to be a circle. We will impose periodic boundary conditions for the circular direction 
and then study what happens when the radius of the circle tends to infinity. A mathematical 
technique, used previously to study the Bose gas in a harmonic oscillator potential 7], and 
a charged Bose gas in a constant mag Mti e Md 0, is usedtoob^ — aptotic 
expansions for a number of thermodynamic properties of the system. A short discussion is 
given at the end of the paper. 

II. THERMODYNAMIC EXPRESSIONS 

We begin with a flat spatial region that is taken to be a finite volume box of dimensions 
L, L 2 , L 3 in the x, y, z directions respectively. The Schrodinger field will be chosen to obey 



2 



periodic boundary conditions on the box walls. We will be interested in the limit where 
L 2 , L 3 — > oo so that the spatial region is topologically S 1 xl 2 . 

Under these conditions the energy levels, with L 2 , L 3 initially finite, are given by 

_ 2n 2 h 2 fn 2 n\ nl\ 
where n, n 2 , n 3 = 0, ±1, ±2, . . .. The total particle number is 

oo oo oo 

N= J2 J2 Yl [e /3 ^"3-M)_i]- 1 (2.2) 

n=— oo n2=— oo n3=— oo 

with /i the chemical potential and (3 = T" 1 the inverse temperature (in units where the 
Boltzmann constant is taken equal to unity). Because the ground state energy vanishes, it 
is clear that 

N 9 r = {e-^-l)~ 1 (2.3) 

gives the ground state particle number. Bose-Einstein condensation viewed as a phase 
transition can only occur if /i = at some temperature, defined as the critical temperature. 
A crucial feature of a finite volume system, or a semi-infinite one such as we have here, is 
that there will be no phase transition characterized by fi = 0. Nevertheless, if we view Bose- 
Einstein condensation as defined by a sudden buildup of particles in the ground state over a 
very small change in temperature, then this can occur even for a finite volume system. An 
important example of where this happens is for the harmonically confined Bose gas 0, E| ■ 
The principle aim of this paper is to study what happens to the system when the size 
of the confining volume is varied. In the present setup, this corresponds to changing the 
circumference of the circle. 

It proves convenient to define two dimensionless parameters x and e by 

. = (2-4) 
pfx = -ex . (2.5) 

If we let L 2 , L 3 become large, then we can approximate the sums over n 2 and n 3 that occur 
in ()2.2j) with integrals and obtain 



N 



^(^)tln[l-e-^]. (2.) 
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The number density is therefore 

N IT 



P 



LL2L3 L 3 x 



In [l - e - x{n2+€) ] . (2.7) 



71= — CO 



We wish to compare our results for finite L with those found in the bulk, or thermo- 
dynamic, limit found when L — > 00. In this limit we can replace the sum over n with an 
integral and find the bulk limit to be the standard result of 

3/2 



Pbulk 



where Li p (e B ) = n p e n6 is the polylogarithm function. The critical value of T at which 

71=1 

[L — in the bulk limit will be defined as To = (3^ 1 where 

3/2 



bid k 



m 



C(3/2) . (2.9) 



Even if we keep L finite, so that the bulk approximation is not made, we can still define 
To as the temperature at which the density p is given by the right hand side of (|2.9jl . We 
will define 



s o = — ttPo, (2.10) 



2n 2 h 2 
mL 2 

and call €q the value of e at the temperature T . Using ()2.9|) and ()2.7|) leads to the relation 



7r 1 / 2 Zo 3/2 « 1/2 C(3/2) = - ln l 1 ~ e- x(n2+e) ] ■ (2.11) 

71= — OO 

The point of doing what we have done is that instead of specifying a particle density, we 
can specify a value for xq instead. The relation (|2.11|) determines e as a function of x for a 
given value of xq. By looking at different values of Xq we can see how results alter as L is 
changed. This is completely equivalent to solving for the chemical potential as a function 
of temperature for a given particle number density, but is better suited to the analytical 
approximations that we will obtain in the next section. 

The result of solving (J2.ll)) numerically for different values of Xq is shown in Fig. ^ It 
can be seen that there is a sharp decrease in e over a very small range of x. From (|2.3|) 
with (|2.4|) this can be seen to correspond to a sudden increase in the number of particles 
in the ground state, and hence to Bose-Einstein condensation. Note that smaller values of 
Xq correspond to larger values of L from (|2.10|) . As x decreases, the value of x at which 




FIG. 1: This shows e plotted as a function of x/xq for different values of xo- 

e undergoes its sudden decrease from values of the order of 1 to values much smaller than 
1, comes closer and closer to the value Xq. In this sense, as L increases we have a situation 
that starts to look more and more like a phase transition. We emphasize that there is no 
true phase transition here because we never attain the value e = for any finite value of L. 

In the bulk limit the specific heat has the familiar peak at the critical temperature To. 
(See Refs. for example. We also show the result in Fig. [21 for comparison.) We 

now wish to see what happens when L is finite. The absence of a phase transition in this 
case indicates that the specific heat should be smooth, unlike the situation for the bulk 
approximation where the derivative of the specific heat is discontinuous at T = T . 

The internal energy is given by 

U= £ E m3 [e /3 ^"3-M)_ 1 ]- 1 . (2.12) 

n,7i2,ri3 

Taking L 2 , L 3 to be large, and replacing the sums over n 2 and ris with integrals again leads 
to 



(V OO 
£ n 2 ln[l-e-^) 
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From this result we can compute the specific heat at constant volume from 



N,V 



After some calculation, it can be shown that 

C v A C 2 



N B BD 



where 



.4 



E {™ 



12 e 



-x(n 2 +e) 



2xn In 



+x 2 n 



x(n 2 +e) 



-1 



B 



C 



D 



- ln l 1 - e~ x{n2+e) 

n=—oo 
oo 

E 



n=— oo 
oo 



xn 



-1 



111 



1 _ g-z(« 2 +e) 



£ e -(n 2 +e)_ 1 



-1 



(2.13) 



(2.14) 



(2.15) 



(2.16) 
(2.17) 

(2.18) 

(2.19) 



Given the numerical solution for e found from (|2.11|) . it is a straightforward task to 
evaluate the sums in ()2.16fJ2~T9"|) and thereby to find the specific heat for given values of x$. 
The results are shown in Fig. for sample values of xq. It can be seen that as xo decreases, 
corresponding to increasing L and becoming closer to the bulk approximation, the specific 
heat maximum shifts closer to x = xq and the shape of the curve near the maximum 
becomes more pointed. This corresponds to what would be expected heuristically, with the 
bulk approximation becoming indistinguishable from the exact result as L becomes larger. 



III. ANALYTICAL EXPANSIONS 

The key ingredient of the method we will use consists of converting the infinite summa- 
tions that occur in the expressions of the previous section to a contour integral using the 
Mellin-Barnes representation for the exponential, 



c+ioo 

/ da T(a)x~ c 

2ni 



(3.1) 
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FIG. 2: This shows Cy /N plotted as a function of x/xq for different values of xq. The result of 
using the bulk approximation is shown for comparison. 

Here c > is a real constant chosen so that the contour lies to the right of all poles in 
any of the resulting expressions. The contour is then closed in the left hand of the complex 
plane and the result evaluated by residues. This method was used previously jlj to obtain 
analytical expressions for the harmonically confined Bose gas. 

We will give a brief description of the method for the sum that occurs in ()2.11)) and then 
simply quote the results for the other sums that are needed. Define 



Start by expanding the logarithm in its Taylor series and then make use of ()3.1|) . This gives 



oo 




(3.2) 



n=— oo 




1=1 



n=—oo 



c+ioo 




(3.3) 



where we have defined 



oo 




(3.4) 



The summation in (|3.4j) converges for dt(a) > 1/2. We therefore choose c > 1/2 in ()3.3|) so 
that all possible poles in the integrand lie to the left of the contour. In the Appendix we 
show that T(a)F(a, e 1//2 ) is an analytic function of a except at a = 1/2 — n for n = 0, 1, 2, . . . 
where simple poles occur with residue ir l l 2 (—e) n /n\. In particular, r(a)i r '(a, e 1//2 ) is analytic 
at a = 0, where ((a + 1) has a simple pole with residue 1. The analytic continuation of 
r(a)F(a,e 1 / 2 ) to a = is given in (|A7[) . We therefore find the asymptotic expansion 

L(e, x) ~ 2 ln[2 sinh(7re 1/2 )] - tt 1/2 ^ C(- - ^)x n ~ 1/2 . (3.5) 

n=0 n ' 

An analysis similar to that that we have just described for L(e, x) may be applied to the 
other sums that occur in (J2.16fJ2~T9*)) . The results are 

A ~ ^ 7r 1 /2 C (5 / 2) a; -i/2_ ^1/2^3/2)^1/2 

+ex {41n[2 sinh(7re 1/2 )] l ' 2 cothfie 1 ' 2 )} 

+T ii/2 £ra c(i/2 - n)a; " +3/2 - (3 ' 6 » 

n=0 

B = -L{e,x) , (3.7) 
C ~ -7r 1/2 C(3/2)a;- 1/2 -7re 1/2 coth(7re 1/2 ) -21n[2sinh(7re 1/2 )] 

+^ 1/2 E^r9c(i/2-n)^/2 ; (3 . 8) 

D ~ vre^^x^ 1 cotl^Tre 1 / 2 ) + tt 1 / 2 ^ C(l/2 - n)x n " 1 ' 2 , (3.9) 

These results should all be viewed as asymptotic expansions in e and x that are valid when 
e and x are both small. It is straightforward to see that this is the case by comparing 
the numerical evaluation of the sums in (|2.1fiti2~T7j) with these analytic expansions. The 
advantage of the method is that when x is small the direct numerical evaluation of the 
sums becomes more difficult because they start to become only slowly converging. When x 
becomes large enough so that the analytical approximations break down (typically x of the 
order of 0.5), the direct evaluation of the sums becomes simple as they converge rapidly. 

It is now straightforward to find an expansion for the specific heat using the analytical 
results of ()3.5H3.9|) in ()2.15|) . We find, after some calculation, 

c ^T(li +aiil/2+a2i+ - ' (3 - 10) 



s 



where the coefficients ot.\ and a.2 are defined by 

ai = 2^ ^| ln[2 sinh (- 1/2 )] - ^ C(3/2)e 1 / 2 tanh^) , (3.11) 
15 C(5/2) 2r • u/ i/2m ^ 

+ — e 1/2 tanh(7re 1/2 ) In [2 sinh(7re 1/2 )] 

+ T C( C^ 1 /^ 6 + ^ C(3/2)C(l/2)etanh 2 ( 7 re 1 / 2 ) . (3.12) 

Only the first few terms are displayed here. This gives us an analytical approximation to 
(12.14)1 that is valid for small x and e. With T fixed, as we let L — > oo we have i -> so 
that only the first term in ()3.10j) survives. This is the familiar value of the specific heat 
maximum from the bulk approximation. If we keep L finite, then for small e it can be shown 
that «i < 0, showing that the peak in the specific heat is smaller than that occurring in the 
bulk limit. This is consistent with what we found numerically as shown in Fig. |21 

Using our analytical approximation ()3.5|) it is possible to obtain an approximate solution 



to (J2.11jl . This will enable us to study the manner in which the chemical potential vanishes 
as L — > oo and the bulk limit is reached. We have called eo the value of e when x = xq 
(corresponding to the temperature T = T defined in Eq. ()2.1()|) ). We then have 

~ 2 In [2 sinh(7re 1/2 )] + tt 1/2 C(1/2) e^ 2 + ■■■ (3.13) 

for small Xq (corresponding to large L). To lowest order in xq we may ignore the second 
term in ()3.13|) and find 



£0 .Wi±^|. (3.14) 



The numerical value of this is eo — 0.0234624, a result that shows up in the numerical 
evaluation of ()2.11|) as xq is decreased. 

______ 1/2 

The next order correction to the result in ()3.14|) will be a term of order x . By expanding 

1/2 

(I3.13|) and working consistently to order x it may be shown that 

-M^K^^M- (3 - i5) 

In a similar way, more terms in the expansion can be evaluated if desired. If we use (I3.15|) 
in (|2.5|) . we can conclude that for large values of L the value of the chemical potential /i at 
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T = T n is, 



A*(r = To) 



2/1 2 



ml 2 



In 2 



1 + V5 



57rmTo 



1/2 



fiC(l/2) 



In 



1 + V5 



+ 



(3.16) 

The next order term inside the braces is of the order (ToL 2 ) -1 . As we let L — » oo, it is seen 
that fi(T = T ) — ► as expected. Our calculation has established in a precise way exactly 
how the bulk limit is reached. 

For temperatures that are close to T = T , a simple extension of the approximate solution 
to (|2.11|) that we have presented gives 



^ln 2 



1 + VS\ | C(3/2)ln( 



l+\/5 
2 



7T^ 



\/5 7T 3 / 2 Xc! 2 



(3.17) 



to lowest order. Again it is straightforward to obtain higher order terms in the expansion if 
needed. 



IV. DISCUSSION AND CONCLUSIONS 

We have studied the ideal Bose gas concentrating on the way that the bulk limit is 
reached. Reliable asymptotic expansions of relevant thermodynamic quantities have been 
obtained and compared with numerical results. The agreement was good in regions where 
the asymptotic expansions are expected to be valid. It was shown how the large volume 
limit of the finite volume system mimics the behaviour of the bulk infinite volume system. 
With more work it is possible to see how the derivative of the specific heat with respect to 
temperature becomes discontinuous as the bulk limit is approached. 

For simplicity of presentation we have only studied the situation for a 3-dimensional 
space with one compact circular direction. There is nothing to stop a similar analysis being 
applied to a space of arbitrary dimension, to more than one compact dimension, or to 
different boundary conditions. It would also be possible to perform a similar analysis for 
the relativistic Bose gas and compare with the results of Shiokawa and Hu [lOj. Finally, it 



would be of interest to study the inclusion o 
possible using the ^-function method in Ref. 



interactions to the analysis. This might be 



111. 
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APPENDIX A: SOME SUM RESULTS 



In this Appendix we outline a derivation of the pole structure of T(a)F(a, e 1//2 ) using the 
same methods as in the main text of the paper. We begin by using the standard integral 
representation [12] for the T-function, assuming initially that dt(a) > 1/2 so that the sum 
in (J3.4|) converges. This gives 

oo 

T(a)F(a,e 1/2 ) = [ dit^e^^t) , (Al) 



o 



where 

oo 

tp{t) = e ~ tn2 • ( A2 ) 

, 13], 



12 



but we do not use this result 



(It can be noted that tp(t) is related to the ^-function 
here.) Because we assume e > 0, and we have <p(t) — > 1 as t — > oo, the integrand of (|A1|) is 
well-behaved for large t no matter what the value of a. Any poles of T(a)F(a, e 1 / 2 ) must 
therefore come from the small t behaviour of the integrand. 

To get the behaviour of the integrand of (|A1J) for small t, we first use fl3.1|) to find that 



oo 

<p{t) = l + 2j2e- tn2 

n=l 

c+ioo 

1 



C— lOO 



1 + — / dz V(z)r z C{2z) . (A3) 



The contour in (|A3|) may be closed noting that there are simple poles at z = 1/2 and z = 
only. Evaluating the residues at these poles results in 

<p(t) ~ vr 1 / 2 t- 1 / 2 (A4) 



121. If we use 



as t — > 0, a result that agrees with what is found using ^-function identities 
VV to mean the pole part of any expression, then it is clear that 

I 

= VV\ vr 1 / 2 jr (a + n- 1/2)- 1 \ . (A5) 

I n=0 U ' 
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(In the last line we have simply expanded e et in its Taylor series and integrated the result.) 



This shows that T(a)F(a, e 1//2 ) has simple poles at a = 1/2 — n with residues tc 1 ^ 2 - — - — for 
n = 0,1,2,.... 

Although this method is a quick way to uncover the pole terms, if we need to know 
T(a)F(a, e 1 / 2 ) at values of a where the result is analytic, then a more elaborate analysis 
must be performed. This has been done by Ghika and Visinescu [14] and in a more general 
case by Ford Ford's result gives us 

oo 

T{a)F{a, e 1/2 ) = vr 1/2 r(a - ^)e 1/2 ~ a + 4r(a) sinvra J dx(x 2 - e)~ a (e 2 ™ - l)" 1 . (A6) 



c 



1/2 



The poles terms found in (|A5|) agree with what is found from this exact result. By letting 
a — > in ()A6|) it is easily seen that 

[T(a)F(a,e^)} \ a=Q = -2 In [2 sinh(vre 1 / 2 )] . (A7) 

We also need the value at a = 1, a result that is easily obtained by noting that 

^ [T(a)F(a, e 1 / 2 )] = -T{a + l)F{a + 1, e 1/2 ) . (A8) 

(This follows from a straightforward differentiation of (|H.4jl .) Making use of (|A7|) and letting 
a — > in ()A8|) gives 

F(l, t l/2 ) = 7te- 1/2 coth(vr e 1/2 ) . (A9) 
This is a standard result j3| and provides a useful check on the results. 
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